Lattice Study on Tensor Glueball in Radiative J /ij) Decay 
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The radiative decay of J/tp into a pure gauge tensor glueball is studied in the quenched lattice 
QCD formalism. With two anisotropic lattices, the mutlipole amplitudes Ei{0), M2(0) and -B3(0) 
are obtained to be 0.114(12) GeV, -0.011(5) GeV and 0.023(8) GeV, respectively. Thus the partial 
decay width T{J/tp — >■ 7G2++) is estimated to be 0.99(22) keV which corresponds to a large branch 
ratio 1.1(2) x 10~^. The phenomenological implication of this result is also discussed. 
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Glueballs are exotic hadron states made up of gluons. 
Their existence is permitted by QCD but has not yet 
been finahy confirmed by experiment. In contrast to the 
scalar glueball, whose possible candidate can be /o( 1370), 
/o(1500) or /o(1710), the experimental evidence for the 
tensor glueball is more obscure. Quenched lattice QCD 
studies predicts the tensor glueball mass to be in the 
range 2.2 — 2.4 GeV [1|-[3j which is also supported by a 
recent 2 + 1 flavor fuU-QCD lattice simulation 4]. In 
this mass region, Mark III Q and BES [a] have ob- 
served a narrow tensor meson ^(2230) (now as /j(2220) 
in PDGQ) in the J/tp radiative decays with a large pro- 
duction rate, whose features favor the interpretation of a 
tensor glueball. However, it was not seen in the inclusive 
7 spectrum [^ by the Crystal Ball collaboration and in 
pp annihilations to pseudoscalar pairs I9l4l6|. So the ex- 
istence of ^(2230)(/j(2220)) needs confirmation by new 
experiments, especially by BESIII experiment with the 
largest J/ip sample. 

It is well known that the production of glueballs is 
favored in J/tp decays because of the gluon-rich environ- 
ment. The radiative decay is of special importance owing 
to its cleaner background. So the production rate of the 
tensor glueball in the decay can be an important crite- 
rion for its identification. The decay has been studied 
only in a few theoretical works |17H19| . In these works 
the tree-level perturbative QCD approach is employed. 
Under certain assumptions the helicity amplitudes of the 
decay are related to the coupling of the two gluons to 
the tensor glueball. This coupling has been determined 
with the quenched lattice QCD ^[20]. Based on results 



of [17l - ll9j the branch ratio is estimated as 2 x 10 ^ [21|, 
but the theoretical uncertainties are not under control. 

In fact, the decay can be investigated directly from the 
numerical lattice QCD studies [23 . l23j . which provides 



first principle calculations from the QCD Lagrangian, es- 
pecially in quenched lattice QCD. Quenched lattice QCD 
can be taken as a theory which only consists of heavy 
quarks and gluons. In this theory amplitudes of the de- 
cay do not have absorptive part because of masses of 
states. Hence, the amplitudes can be directly calculated 
in the theory in Euclidian space. It should be noted that 
it is still a challenging task for the full-QCD lattice study 
of the decay because glueballs can be mixed with states 
of light quark-pairs. Nevertheless, the study of the decay 
in quenched QCD will give important information about 
nonperturbative properties of glueballs. 

At the lowest order of QED, the amplitude for the 
radiative decay J/tp -^ ^02++ is given by 
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where q = Pi ~Pf is the momentum of the real photon, r, 
r^ and vq are the quantum numbers of the polarizations 
of J/tp, the photon and the tensor glueball, respectively. 
ep((f, r^) is the polarization vector of the photon and j^ 
is the electromagnetic current operator. The hadronic 
matrix element appearing in the above equation can be 
obtained directly from lattice QCD calculation of corre- 
sponding three-point functions. On the other hand, these 
matrix elements can be expressed (in Minkowski space- 
time) in terms of multipole form factors as follows: 



+a^Es{Q^) + <Ci(Q2) + a^C2(Q2) 
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where af are Lorentz-covariant kinematic functions of 
Pi and pf (and specific polarizations of the states) , 
whose explicit expressions can be derived exactly |24l.l25j. 
EiiQ^), AUQ^), EsiQ^), CiiQ^) and C2{Q^) are the 
form factors which depend only on Q^ = —{Pi ~Pf)^- 
Since Ci(Q^) and C2{Q^) vanish at Q^ = 0, we focus on 
the extraction of the first three which are involved in the 



calculation of the decay width T{J/tp — > 7G2++) as 



r = 
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(|i?i(0)p + |M2(0)|2 + |i?3(0)H, (3) 



where a is the fine structure constant, pj is the photon 
momentum with Ip^,] = (Mj^^ - Mg)/(2Mj/^). 

We use the tadpole-improved gauge action [l| to gen- 
erate gauge configurations on anisotropic lattices with 
the aspect ratio ^ — as/ at — 5, where as and at are the 
spatial and temporal lattice spacing, respectively. Two 
lattices L^xT^S^x 96(/3 = 2.4) and 12^ x 144(/3 = 2.8) 
are applied to check effect of the finite lattice spacings. 
The relevant input parameters are listed in Table|Tl where 
as values are determined from r^^ = 410(20) MeV. Since 
glueball relevant study needs a quite large statistics, the 
spatial extensions of both lattices are properly chosen 
to be '^ 1.7 fm according to the study of finite volume 
efl^ect study (FVE) of Ref. ^], which is a compromise 
of the computational resource requirement and negligi- 
ble FVEs both for glueballs [3| and charmonia. In the 
practice, we generated 5000 configurations for each lat- 
tice. The charm quark propagators are calculated using 
the tadpole- improved clover action for anisotropic lat- 
tices [23, [231 with the bare charm quark masses set by 
the physical mass of J/ip, Mj/^ — 3.097 GeV, through 
which the spectrum of the 15' and IP charmonia are well 
reproduced [2J|. In practice, disconnected diagrams due 
to the charm and quark- antiquark annihilation are ex- 
pected to be unimportant according to the OZI rule and 
therefore are neglected in the calculation of relevant two- 
point and three-point functions. 



TABLE I: The input parameters for the calculation. Values 
for the coupling /3, anisotropy ^, the lattice spacing as, lattice 
size, and the number of measurements are listed. 



/? 


C 


as (fm) 


Las (fm) 


L^ xT 


Accent 


2.4 
2.8 


5 
5 


0.222(2) 
0.138(1) 


1.78 
1.66 


8'' x96 
12=* X 144 


5000 
5000 



The calculations in this Letter are performed in the 
rest frame of the tensor glueball. One of the key issues in 
our calculation is to construct optimal interpolating field 
operators which couple dominantly to the pure gauge 
tensor glueball. This is realized by applying completely 
the same scheme as that in the calculations of glueball 
spectrum [2, Sj- On the cubic lattice, a tensor (J = 2) 
state corresponds to the E and T2 irreducible represen- 
tations of the lattice symmetry group O. So we build 
the E and T2 operators from a set of prototype Wilson 
loops. By using different gauge-link smearing techniques. 



an operator set {(/)«,« = 1,2, ...,24} of 24 different 
gluonic operators is constructed for each component of 
the T2^ and i?++ representations, where the superscript 
i labels the three components of T2 and two compents 
of E. Finally for each component, an optimal operator 
<!>(*)(/;) = Y^ Va't'a {t) for the ground state tensor glueball 
is obtained with the combinational coefficients Va deter- 
mined by solving the generalized eigenvalue problem 



C'^'Hiz))v(-^) = e-*°™(*")C'(*)(0)v(-«\ 



(4) 



at to = 1, where C^^^{t) is the correlation matrix of the 
operator set, 

T 

In addition, the glueball two-point functions are normal- 
ized as 

C\t) = i^($W(t + r)<i>Wt(r)) 
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where \Ti) refers to z-th component of the Tj"*"*" and E^^ 
glueball state. We are assured that C*(t) can be well de- 
scribed by a single exponential, C{t) — We^^^'^*, with 
W usually deviating from one by few percents. It should 
be noted that, due to the SO{3) rotational symmetry is 
broken on the lattice with a finite lattice spacing, and 
consequently the masses of T2 and E glueballs are not 
necessarily the same, even though they converge to the 
same tensor glueball mass in the continuum limit when 
the rotational invariance is restored. However, with the 
two lattice spacings we used in this Letter, we observe 
that the difference of the two masses are not distinguish- 
able within errors, which implies that the effects of the 
rotational symmetry breaking are not important. So in 
the following, we neglect this symmetry breaking and as- 
sume that the five components of the T2 and E and that 
of the corresponding spin-two state can be connected by 
a normal transformation. 

We calculate the three-point functions in the rest frame 
of the tensor glueball with J/tp moving with a definite 
momentum pj = 2Trn/Las, where ft ranges from (0, 0, 0) 
to (2, 2, 2). In order to increase the statistics additionally, 
for each configuration we calculate T charm quark propa- 
gators Spix, i; 0, r) by setting a point source on each time 
slice r, which permits us to average over the temporal di- 
rection when calculating the three-point functions 



T-1 



rgL('f;i/,0 - ^ 5] 5]e-^«~-^>«(i/ +r)J^(y,i + r)Ov,,(0,r)) 



E 



T=0 a 

g_MT(t/-t)g-£v(f)i 



2MTV32Ev{q) 



{0\<i>^^\0)\T.)m\M0)\V{q,r)){V{q,r)\OU0)\0), 



(7) 



;=•■■. 


000 ' — • — ' 
100 ^^ 

110 — -- 

111 -^ 

200 -^ 

210 ^— ^ 

211 ^--^ 

220 -^ 

221 ^— ^ 

222 ^--^ ■ 



FIG. 1: The effective energy plot for J/tp with different spa- 
tial momenta. From top to bottom are the plateaus for 
momentum modes, p — 2Tin/L, with ft = (2,2,2), (2,2,1), 
(2,2,0), (2,1,1), (2,1,0), (2,0,0), (1,1,1), (1,1,0), (1,0,0), 
and (0,0,0). 



where J^{x) — c.{x)^^c(x) is the vector current operator, 
Ov,j = c'^jC the conventional interpolation field for J/V', 
and the summation in the last equality is over all the 
possible states with different polarizations. In the rest 
frame of the tensor glueball, the momentum of the initial 
J /ip is the same as that of the current operator, say, 
pv = q. The vector current J^{x), which is conserved 
in the continuum limit, is no longer conserved on the 
lattice and requires a multiplicative renormalization. The 
renormalization constant of spatial components of the 
vector current is determined to be Zy = 1.39(2) for 



13 = 2.4 and Z'y' = l-ll(l) for /3 = 2.8 [22| using a 
nonperturbative scheme |23| . 

The matrix elements (Ti| J^(0)|y(g, r)) can extracted 
from the above three-point functions along with the two- 
point function of the glueball C*(i) and that of J/ij} 

Tf\q.t), 



Tfm = Y^e 



-''"-^(0|Ov,,(f,t)O^,,(0,0)|0), (8) 



which provide the information of Mt, Ev{q) and the 
other two matrix elements. According to Eq. ^ one has 
approximately 



(0|$«(0)|l^,) w V2Mt1^3. 



(9) 



Mt and Ey{q) can be determined precisely from the 
two-point functions. Figure [T] shows the nice effective 



energy plateaus of Jf-ijj for typical momentum modes 
at /3 = 2.4. We also check the dispersion relation of 
J/f/) and find the largest deviation of squared speed of 
light c^ from one is less than 4%. The matrix elements 
(V{q,r)\Oy A0)\0} are included implicitly in the three- 
point and two-point functions and can be cancelled out 
by takeing a ratio Ri^fj_j{q,t) for some specific {i,fi,j} 
combinations. 
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which should be insensitive to the variation of i in a time 
window. As such the desired matrix elements can be 
derived by the fit form 

i?,,^,,(q,i) = J2iT^\J^mVi<l,r))eJiq,r)+Sf{t), (11) 



where e{q, r) is the polarization vector of J/^p and 5f(t) — 
ae"™* accounts for the residual contamination from ex- 
cited states. 

In the data analysis, the 5000 configurations are di- 
vided into 100 bins and the average of 50 measurements 
in each bin is taken as an independent measurement. 
For the resultant 100 measurements, the one-eliminating 
jackknife method is used to perform the fit for the matrix 
elements (Mt and Ey determined from two-point func- 
tions are used as known parameters). Generally speak- 
ing, the time separation t and tj —t should be kept large 
for the saturation of the ground state, but we have to fix 
tf —t ^ 1 because of the rapidly damping of the glueball 
signal. Fortunately this is justified to some extent by 
the optimal glueball operators which couple almost ex- 
clusively to the ground state. The second step of the data 
analysis is to extract the form factors Ei{Q'^), M2{Q'^) 
and EsiQ^) at different Q^ = 2Ev{q)MT - M^ - M| 
according to Eq. ([2]). Since the matrix elements are mea- 
sured from the same configuration ensemble, we carry out 
a correlated data fitting to get these three form factors 
simultaneously with covariance matrix constructed from 
the jackknife ensemble described above. The symmetric 
combinations of the indices {i,/^, j} and the momentum 
q which gives the same Q^ are averaged to increase the 
statistics. In order to get the form factor at Q^ = 0, we 




FIG. 2: The extracted form factors -Ei(Q^) M2{Q^) and 
E'i{Q^) in the physical units. The upper panel is for /3 — 2.4 
and the lower one for /? = 2.8. The curves with error bands 
show the polynomial fit with F^{Q'^) = Fi(0) + a^Q^ + hQ'^ . 



carry out a correlated polynomial fit to the three form 
factors from Q^ 



-0.5 GeV^ to2.7GeV^ 



^.(Q')=i^»(0) + a,;g2 + 6,g4 



(12) 



where Fi refers to Ei , M2 or E-^ . Figure [2] shows the 
final results of Fi{Q-^) for /3 = 2.4 (upper panel) and 
/3 = 2.8 (lower panel), where the data points are the 
calculated values with jackknife errors, and the curves 
are the polynomial fits with jackknife error bands. The 
concrete numerical results of the tensor glueball mass Mq 
and the three form factors on the two lattice are also 
shown in TABLE. Illl Note that renormalization constant 

(s) — 

Zy of the spatial components of the vector current is 
applied to the final numerical values. 

The last step is the continuum extrapolation using the 
two lattice systems. By performing a linear extrapola- 
tion in Cj , the continuum limit of the three form factors 
are determined to be £^1(0) = 0.114(12) GeV, M2(0) = 
-0.011(5) GeV and ^3(0) = 0.023(8) GeV, respectively. 
Note that there is a pattern \Ei{0)\ > |M2(0)| ~ |^3(0J|, 
hence the decay width is dominated by the value of 
-Ei(O). For the continuum value of the tensor glueball 
mass, we take Ma = 2.39(12) GeV from Ref. [3]. Thus, 
according to Eq. ([3]), we finally get the decay width 
r{J/ip -^ -1G2++) = 0.99(22) keV. Using the reported to- 
tal width of J/^, Ftot = 92.9(2.8) keV, the corresponding 
branching ratio is 



r(j/^ ^ 7G2++)/rtot - 1.1(2) X 10- 



(13) 



2.4 2.360(20) 


0.142(07) 


-0.012(2) 


2.8 2.367(25) 


0.125(10) 


-0.011(4) 


cx) 2.39(12) f3] 


0.114(12) 


-0.011(5) 



TABLE II: Listed in the table are the tensor glueball masses 
Mt and the form factors Ei, M2 and -E3, calculated on the 
two lattices with /3 = 2.4 and 2.8, respectively. Also shown are 
the continuum extrapolation and the resultant partial width 
T. 

P MT(GeV) El (GeV) Ma (GeV) E3 (GeV) r(kcV) 

0.012(2) 1.46(18) 

0.019(6) 1.17(20) 

0.023(8) 0.99(22) 



The determined branching ratio is rather large. We 
admit that the calculation is carried out in the quenched 
approximation, whose systematical uncertainty cannot 
be estimated easily without unquenched calculations. A 
recent fuU-QCD lattice study of the mass spectrum of 
glueballs in Q indicates that there is no substantial cor- 
rection of the masses of the scalar and tensor glueball. 
Based on this fact one can expect that the unquenching 
effects may not drastically change our result. 

In experiments, the narrow state /j(2220) observed 
by Mark III and BES in the J /ijj decay was once inter- 
preted as a candidate for the tensor glueball. But the 
analysis of the processes pp -^ '!T'k{KK) yields no indica- 
tion of the narrow /j(2220) and sets an upper bound for 
the branch ratios Br{fj -^ pp)Br{fj -^ Tnr,KK) (see 
the review article Ref. |23| and the references therein). 
Combining this with the results of Mark III and BES, 
a lower bound for the branching ration_is obtained to 
be Br(J/V' ^ 7/j(2220) > 2.5 x lO'^ fl], which seems 
compatible with our result. However, BESII with sub- 
stantially more statistics does not find the evidence of a 
narrow structure around 2.2 GeV of tttt invariant mass 
spectrum in the processes J/-0 — >■ 77777 ^9]. Recently, 
based on 225 million J/ip events, the BESIII Collabora- 
tion performs a partial wave analysis of J/ip — >■ 77777 and 
also finds no evident narrow peak for /j(2220) in the rji] 
mass spectrum [30]. So the existence of fj {2220) is still 
very weak. Another possibility also exists that the ten- 
sor glueball is a broad resonance and readily decays to 
light hadrons. Our result motivates a serious joint analy- 
sis of the radiative J/^p decay into tensor objects in VV, 
PP, pp and 47r final states (where V and P stand for vec- 
tor and pseudoscalar mesons, respectively), among which 
VV channels may be of special importance since they are 
kinematically favored in the decay of a tensor meson. 

To summarize, we have carried out the first lattice 
study on the Ei, M2 and E3 multipole amplitudes for 
J/ip radiatively decaying into the pure gauge tensor glue- 
ball G2++ in the quenched approximation. With two 
different lattice spacings, the amplitudes are extrapo- 
lated to their continuum limits. The partial decay width 
and branch ratio for J/ip -> 7G2++ are predicted to be 
0.99(22) keV and 1.1(2) x lO'^, respectively, which imply 
that the tensor glueball can be copiously produced in the 
J/ip radiative decays if it does exist. To date the exis- 
tence of /j(2220) needs confirmation and a broad tensor 



glueball is also possible. Hopefully, the BESIII data will 
be able to clarify the situation. 
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